Non-Universality in Random Matrix Ensembles with Soft Level Confinement 
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Two families of strongly non-Gaussian random matrix ensembles (RME) are considered. They 
are statistically equivalent to a one-dimensional plasma of particles interacting logarithmically and 
confined by the potential that has the long-range behavior V(e) ~ |e| Q (0 < a < 1), or V(e) ~ 
In 2 |e|. The direct Monte Carlo simulations on the effective plasma model shows that the spacing 
distribution function (SDF) in such RME can deviate from that of the classical Gaussian ensembles. 
For power-law potentials, this deviation is seen only near the origin e ~ 0, while for the double- 
logarithmic potential the SDF shows the cross-over from the Wigner-Dyson to Poisson behavior in 
the bulk of the spectrum. 
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The classical theory of Random Matrices (RMT) devel- 
oped by Wigner, Dyson and Mehta [Q provides a statis- 
tical description of the energy levels in a variety of quan- 
tum chaotic systems. In this way, one of the simplest 
statistical characteristics is the probability distribution 
P(s) of the spacing between nearest-neighbor eigenval- 
ues. In the framework of the classical RMT, the spacing 
distribution function (SDF) follows very closely a univer- 
sal curve known as Wigner surmise [lj. Its most impor- 
tant characteristic is the vanishing of P(s) at s = 0, that 
demonstrates level repulsion. 

In contrast, for classically non-chaotic systems, the 
random energy levels are described by another univer- 
sal distribution, the Poisson statistics, which assumes all 
levels to be uncorrelated. 

Both universal statistics are realized in a disordered 
system of non-interacting electrons. The metal phase 
that exists for relatively weak disorder, was proved ||] 
to be described by the Wigner-Dyson statistics while the 
level statistics in the insulator phase is close to the Pois- 
son distribution. 

The transition between this two phases known as the 
Anderson transition, has much in common with the criti- 
cal phenomena in second order phase transitions and can 
be described by the scaling approach ||]. Using scaling 
arguments, one can show that in the critical region near 
the Anderson transition, there should exist a third uni- 
versal statistics [Q. The detailed scaling analysis done 
recently || showed this statistics to be drastically differ- 
ent from both the Wigner-Dyson and the Poisson statis- 
tics, the corresponding spectral correlation functions be- 
ing characterized by nontrivial exponents related to the 
correlation length exponent v. 

It is of great interest to see if a description of the criti- 
cal statistics in terms of random matrices is still possible. 
Clearly, if this is the case, the corresponding matrix en- 
sembles must be of a completely different kind from the 
ones belonging to the Wigner-Dyson universality class 



or from those leading to the Poisson statistics. Our in- 
terest in this direction was prompted by another recent 
discovery || of a new family of random matrices, a one- 
parameter solvable model, that displays a cross-over in 
the spacing distribution from a highly correlated Wigner- 
Dyson to a completely uncorrelated Poisson distribution 
when the parameter is varied. 

In this paper we have tried to establish if the non- 
classical behavior of level correlations is a generic feature 
shared a by broader class of random matrix models. We 
will consider two strongly non- Gaussian ensembles of ran- 
dom matrices. We will show that the first one breaks the 
Wigner-Dyson universality only locally, in the center of 
the spectrum. It will nevertheless allow us to understand 
better the second class of models which is similar to the 
exactly solvable model studied in Ref. Q. Here we will 
show that this class of random matrices indeed breaks 
the Wigner-Dyson universality globally and displays a 
cross-over to a Poisson-like distribution. 

Let us consider a physical system described by N x N 
random matrix H whose eigenvalues {e n },n = l,...N 
will be also randomly distributed. Within the maximum 
entropy ansatz for describing the eigenvalue distri- 
bution at a given mean level density p(e), we can use the 
effective plasma model introduced by Dyson [jjj, where 
the joint probability density function P({e}) is mapped 
onto the Gibbs distribution of a classical one-dimensional 
plasma of fictitious particles with a pair-wise logarithmic 
repulsion — In |e„ — e m \ and a one-particle potential V(e) 
to keep the system confined: 



P({e}) = Z- 1 exp[-PH({e n })}, 
H({e n }) = -5>| ei - ej -|+£y( f 
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Here Z is the partition function and 0, that plays the role 
of an inverse temperature in the corresponding Gibbs en- 
semble, is related to the symmetry of the original random 
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matrix ensemble and is equal to 1,2 or 4 for orthogonal, 
unitary and symplectic ensembles, respectively |Q. 

When the confining potential is quadratic, the model 
is the Gaussian ensemble studied by Wigner, Dyson and 
Mehta and the corresponding SDF is very close to the 
Wigner surmise. The inclusion of higher powers of e 2 
was shown || to make no effect on the SDF in the limit 
N — ► oo. Until very recently the confining potential was 
believed to be irrelevant for level correlations in the ther- 
modynamic limit N — > oo. 

However, a recent work || has demonstrated that it is 
not the case. For some specific one-parameter family of 
confining potentials V p (e), the exact solution in terms of 
the non-classical q-orthogonal polynomials was found. It 
showed the level correlations to deviate from the conven- 
tional Wigner-Dyson form as the parameter p increases, 
the SDF approaching the Poisson distribution for large 
values of p. This was associated with the asymptotics of 
the confining potential V p (e) ~ In 2 |e| for |e| 3> 1 that 
is an extremely "soft" confinement as compared to the 
Gaussian confining potential V(e) — e 2 . 

In this paper two main questions will be addressed. 
The first one is how soft should the confining potential 
be in order to see deviations from the classical Wigner 
surmise. The second question is whether the exact solu- 
tion found in Ref. || represents the generic features of all 
models with the double-logarithmic long-range behavior 
of the confining potential. 

In order to answer the first of these questions, we con- 
sider a family of power-law potentials: 



V(e) 



A, 
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where A > and a > are two constant parameters. For 
a = 2 Eq.(||) reduces to the Gaussian quadratic confine- 
ment, while in the limit a — > the combination of such 
power-law potentials reproduces the double-logarithmic 
potential V(e) = In 2 |e| = lim Q ^ [a- 2 (|e| Q - l) 2 ]. 

We will see that the level statistics exhibits two sharp 
transitions when the parameter a decreases. The first one 
occurs at a = 1 and it is connected with the break-down 
of translational invariance in the eigenvalue space that is 
present for a > 1 in the limit N — > oo. For a < 1, the 
SDF shows a non-classical, a-dependent behavior only 
near the center of the spectrum e = 0. The second critical 
value is a = 0. For confining potentials that increase 
only logarithmically, the SDF turns out to deviate from 
the Wigner-Dyson form everywhere in the bulk of the 
spectrum. 

The first transition is seen already within the mean- 
field approximation suggested by Dyson pj|. Let us de- 
fine p(e) = y^i S(e — Ci). By substituting this definition 
into the Eq. (Q) one obtains the continuous version of the 
energy functional Ti\p] in terms of p(e). The extremum 
of this functional corresponds to an equilibrium of the 
effective plasma expressed by the equation: 



de'(p(e')) \n\e-e'\ = V(e) + c, 



(4) 



where (p(e)) is the mean density, and the Lagrange mul- 
tiple c is to be found from the normalization condition 
f(p(e))de = N. 

Such a mean-field (MF) approximation completely dis- 
regards the entropy part of the free-energy functional and 
is exactly applicable only for (3 — oo. However, the long- 
range nature of the pair- wise interaction in Eq. (^) makes 
the MF approximation applicable in the bulk of the spec- 
trum even at finite /?, since for the class of confining po- 
tentials of Eq.(|j|) the relative entropy contribution van- 
ishes in the thermodynamic limit as (iVlnTV) -1 . 

The solution /9mf(c) to the MF Eq.(||), confined to the 
region —D < e < D, can be found using the Cauchy 
method fl0|| and is given by: 
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where e + = e + iO and the band-edge D is to be found 
from the normalization condition. 

For a > 1 (strong confinement) the main contribution 
to the integral in Eq.(||) is made by the region £ ~ D. 
In the thermodynamic limit N — > oo the band-edge is 
also divergent D — > oo. Therefore, for any fixed |e| <ti 
D, one can neglect the e-dependence in the integrand of 
Eq.(||). Then the mean level density tends to a constant 
p ~ TV 1 ^ 1 /", signaling the translational invariance in the 
e-space. 

However, for a < 1 (weak confinement), the integral in 
Eq.(||) is convergent even in the limit D — * oo. The corre- 
sponding limiting function Pmf( £ ) k N"" 1 can ^ e easily 
found as the limit z = e/D — > Oof the exact solution 
Pmf(c) to Eq.(§): 



Pmf(c) = AC a 



Here C a 
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i>/2)r(i+a/2) > F(a,b;c;x) is a hypergeo- 
metric function, and the band-edge is given by D = 

2 ( NV 2 ( a /2) \ 1/a 

z \ 2Ar( a ) ) ■ 

Thus for a < 1 the mean density Eq.(||) shows the 
lack of translational invariance in the large- N limit and 
is singular at e = 0. 

This singularity, however, appears only in the MF ap- 
proximation. An exact treatment for /3 = 2 that is based 
on the representation in terms of orthogonal polynomials, 
prf shows the value (p(0)) to be finite: 



<P(0)> = 



(A/n 



l/a 



(2/a)r(l/a) ^ 



r(i + l/2) 



r(i + i) 



2/a 
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Thus in case of weak confinement the MF approxi- 
mation fails to describe the mean level density near the 
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origin. It is natural to suppose that all the level correla- 
tion functions will also have a non-classical form in this 
region. 

In order to study the correlation functions and in 
particular the SDF, we have exploited the Coulomb 
plasma analogy and carried out systematic Monte Carlo 
(MC) simulations on the one-dimensional classical sys- 
tem whose probability distribution is given by Eqs.(|l],||). 
As a check that this method works and is numerically 
accurate we have first studied the three Gaussian ensem- 
bles whose density, two-point correlation functions and 
spacing distribution are exactly known 0]. For these 
systems MC turned out to work extremely well for each 
of these quantities. For the power-law potential Eq.(||) 
for a < 1, we have carried out simulations of systems 
up to N — 200 particles. The simulations are very sta- 
ble even for smaller N and we have typically worked with 
N = 100. The evaluation of the mean density is straight- 
forward. In Fig. |l| we plot this quantity for a — 0.5 and 
(3 = 1, 2, 4. The Monte Carlo result agrees very well with 
Pmf found from Eq.(0) except around the origin, where 
the simulation is more accurate and correctly gives a fi- 
nite density at e = 0. For (3 = 2 the Monte Carlo value 
coincides with that found from Eq.(Q). The simulations 
with different numbers of particles illustrate another im- 
portant property of the particle density for weak con- 
finement (a < 1), that is the "incompressibility" of the 
core of the particle-density distribution. In contrast to 
the a > 1 case, for a < 1 the confining potential is so 
weak that it does not "compress" particles in the core 
region near the origin. On adding more particles to the 
system, these get positioned about the wings of the distri- 
bution, rather than distribute themselves homogeneously 
throughout the spectrum, as in the case of strong con- 
finement (a > 1). The particle density p(e) in the core 
region is almost independent of the number of particles 
but depends on the inverse temperature (3. 

The latter dependence is also a characteristic feature 
of the weak confinement. For strong confinement, the 
/3-dependence is present only in 1/N corrections to the 
mean density and thus is negligible. It leads, in particu- 
lar, to the independence of the mean level density of the 
symmetry of the Hamiltonian. For random matrix en- 
sembles with weak confinement considered here, all the 
/?-dependence is "accumulated" in the core region near 
the origin that contains a few levels on the average. 

The MC evaluation of the SDF is, in principle, also 
straightforward. However, in order to compare it with 
the Wigner surmise we need to rescale the particle posi- 
tions e so that the average spacing between two adjacent 
ones is one. This is known as an "unfolding procedure" 
and is always used in numerical calculations of spectral 
correlations [l2|] . It consists in introducing the new vari- 
able a instead of e according to a map: 



a(e)= / W))de'. 
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The mean density is trivially unity as a function of this 
variable. 

In order to study the SDF in the bulk, we use the 
MF solution Eq.(||) for unfolding according to Eq.(||). 
We checked that the obtained unfolded mean density is 
consistently equal to one, except close to the origin and 
the band-edge. The unfolded spacing P(cr) turned out, 
within our numerical accuracy, identical to the Wigner 
surmise for any a. Therefore, in the bulk of the spectrum, 
the Wigner-Dyson universality holds for the power-law 
weakly confining potentials. 

However this universality is broken around the origin. 
To show this we consider a reference particle fixed at 
the origin. The unfolded spacing must be evaluated in a 
different way here, since the MF density is not accurate. 
Therefore we perform the unfolding by computing 



P(a) 



P(e) 



Pie) 
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e(a) 



where the function e(a) is obtained by inverting numeri- 
cally Eq. (||) and using for p(a) the density evaluated by 
MC simulations. The result is shown in Fig. |^, where we 
plot the unfolded SDF for few values of a < 1 and (3=1 
in comparison with the classical spacing of the Gaussian 
ensemble. We can clearly see that P(cr) for small a does 
not follow the Wigner-Dyson universal behavior a 13 and 
starts out roughly like ct' 9 /". If we assume that the new 
variable a is proportional to e a , as obtained from Eq.(||) 
using p(e) oc e a ~ l , this behavior would correspond to 
P(e)/p(e) (x e 13 . Notice also that the decay of the SDF 
for s > 1 depends on a and is slower than that for the 
Wigner-Dyson distribution. We conclude that for the 
power-law weak confining potential, the Wigner-Dyson 
universality is broken only locally around e ~ 0. This 
conclusion is also reached for (3 = 2, using the indepen- 
dent method of orthogonal polynomials. jjlj 

Now we consider the second class of random matrices, 
with the confining potential that behaves asymptotically 
like V(e) oc In 2 |e|. Since our goal is to study the eigen- 
value correlations in the bulk of the spectrum, we choose 
for numerical simulations the regularized confining po- 
tential that is equal to zero at the origin: 



V(e) 



A 



ln 2 (l + B|e| 



(10) 



where A and B are parameters of order 1. 

In Fig. | we show the bulk SDF for A = 1, 0.5, 0.2, 0.1, 
for the orthogonal symmetry ((3 = 1), together with 
the spacing distribution of the corresponding Gaussian 
ensemble for comparison. We can clearly see that for 
small enough A the spacing distribution departs from 
the Wigner distribution and shows an incipient tendency 
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to become more Poisson-like when A is further reduced. 
Similar deviations from the Gaussian ensemble occur also 
for the unitary and symplectic case ((3 = 2,4). 

This is very similar to the cross-over found analyti- 
cally in Ref. for the exactly solvable model with the 
double-logarithmic long-range behavior of the confining 
potential. We can conclude, therefore, that the cross-over 
is indeed not an exclusive property of the exactly solv- 
able model and is more likely a generic feature shared 
by all the random matrix ensembles with the double- 
logarithmic asymptotics of the confining potential. 

The cross-over in the spacing distribution displayed by 
this family of random matrices is remarkably similar to 
the transition observed in exact numerical calculations 
P,|6|Jl^1 on finite disordered systems going through the 
Anderson transition. 



FIG. 2. Nearest-neighbor spacing distribution in the mid- 
dle of the spectrum for /? = 1 and different values of a. The 
a = 2 case corresponds to the Gaussian Orthogonal Ensem- 
ble. 

FIG. 3. Nearest-neighbor spacing distribution for the loga- 
rithmic confining potential, measured in the bulk of the spec- 
trum at P — 1. By decreasing the parameter A, the spacing 
deviates from the universal Wigner-Dyson distribution ap- 
proaching the Poisson distribution (both also plotted). 
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FIG. 1. Particle density for power-law potential with 
a = 0.5: The Monte Carlo results for /3 = 1,2,4 are plot- 
ted vs. the MF density. 
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